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^ ^' Abstract. Log-periodic amplitudes of the surface magnetization are calculated 

. analytically for two Ising quantum chains with aperiodic modulations of the 

couplings. The oscillating behaviour is linked to the discrete scale invariance of the 
perturbations. For the Fredholm sequence, the aperiodic modulation is marginal and 
the amplitudes are obtained as functions of the deviation from the critical point. For 
^ . the other sequence, the perturbation is relevant and the critical surface magnetization 

OO ' is studied. 

m 
o 

O ' 1. Introduction 



in 



The possible occurrence of log-periodic critical amplitudes in sytems with a discrete 
scale invariance has been known since the early days of the renormalization group in 
statistical mechanics [1-3]. 
I Under length rescaling by a factor of 6, in the case of a one-parameter 

Q ' renormalization group, the linearized recursion relation for the singular part of some 
O ■ observable F{9) takes the form: 



X : i^sin.(^) = - Kin^if^O) , (1.1) 

where a and fi depend on the dilatation factor b. Looking for a power-law solution 
under the form -Fsing(^) = ^(^) obtains x = Ina/ln^u and A{9) = A{fj,9). 

The condition on the amplitude is satisfied by a periodic function of In 9 with period 
In/x. When the scale invariance is continuous, b remains arbitrary and cannot enter 
the amplitude which then reduces to a constant. With a discrete scale invariance on 
the contrary, a dependence on the fixed scaling factor is allowed and the amplitude 
may be log-periodic. Such a behaviour may be ascribed to complex exponents since a 
Fourier component may be written as Re{^^+'^ } = cos(x'ln^) 9^. 

Oscillating amplitudes indeed appear in systems with built-in discrete scale 
invariance. One may mention low-dimensional dynamical models for the transition 
to chaos [4-8] , wave propagation on discrete fractals [9] , superconductive transition on 
Sierpinsky networks [10], spin systems on hierarchical lattices [11]. 

More recently, the same behaviour have been also observed in systems like fracture 
in heterogeneous media [12], foreshock activity preceding major earthquakes [13] and 
diffusion-limited aggregation [14], where the discrete scale invariance is not apparent. 
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Such log-periodic modulations might be of great practical importance in refining the 
prediction of earthquakes [13,15]. 

The purpose of the present work is to study analytically such oscillating amplitudes 
which appear in the surface critical behaviour of aperiodic Ising quantum chains with 
Hamiltonian 

^ oo 

^ = -P E + ^k^>l+i] ' ^fe = -^^'^ ' (1-2) 
^ fc=i 

where the a^^s are Pauli spin operators, A is the unperturbed coupling and r the 
modulation amplitude of the couplings. The aperiodic sequence, /jt, = or 1, is 
generated via an inflation rule, leading to a discrete scale invariance for the modulation. 

Previous work on spin systems involved a linear approximation for the recursion 
relation (1.1), leading to errors of several decades in the oscillation amplitudes [11]. 
Here the recursion can be treated without any approximation and we obtain exact 
results for the log-periodic amplitude of the surface magnetization with two different 
modulations. We first complete the study of a marginal sequence for which the leading 
critical behaviour is known [16]. Then we consider a relevant aperiodic modulation, 
leading to a nonvanishing surface magnetization at the bulk critical point. 



2. Marginal surface magnetization: the generalized Fredholm sequence 

The generalized Fredholm sequence, which is the characteristic sequence of the powers 
of m, follows from substitutions on the letters A, B and C: 

A SiA) = ABCC ■■■ C 
B^S{B) ^BC C C ■■■ C 
C S{C) = C C C C ■■■ (7 

m 

With words of length m = 2, one recovers the usual Fredholm substitution [17]. 
Starting with A and associating = to A and C and = 1 to S, for m = 2, 
one obtains the sequence: 

OIIOIOOOIOOOOOOOIOOO--- (2.2) 

More generally, f\. = l when k = mP + 1 so that — X]i=i fk satisfies the recursion 
relation n^i_^^ — ni_^_^ + 1 {q = l,m; = 0) which can be iterated to give n^^ ~ InL. 
The density of perturbed couplings, Uj^^/L, asymptotically vanishes: the Fredholm 
aperiodicity is a surface extended perturbation which does not affect the bulk critical 
properties of the system [16]. The critical coupling, in particular, remains at = 1. 

An heuristic scaling argument has been recently developped by Luck [18], 
extending to aperiodic perturbations the Harris criterion for random systems [19]. 
According to Luck's criterion, the relevance of the perturbation is governed by the 
crossover exponent (p = 1 + v{uj — 1) which involves the correlation length exponent v 
(equal to 1 for the = 1 + 1 Ising model) and the wandering exponent uj of the 
aperiodic sequence [20, 21]. For the Fredholm sequence, a; = so that = and the 



(2.1) 
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perturbation is marginal: the surface exponents vary continuously with the modulation 
amplitude r [16]. 

The surface magnetization may be expressed as a series involving products of 
the couplings A^^, [22] 



oo J 



™.=(i+Eri^r)"'' (2-3) 

j=i fc=i 

which may be rewritten as: 



rris 



[S{X, r)] , S{X, r)^Yl ^'''r-'^' , n, ^ . (2.4) 



With T(A, r) = S'(A, r) — 1 — A ^, the following recursion is obtained [16]: 

(A^ - 1) T(A, r) = r-2(A-2 - A'^'") + r'^iX^"^ - 1) T(A™, r) . (2.5) 

In the recursion process, A is changed into R[X] = A"^. Defining 9 = InA^, one obtains 
a linear renormalization in the new variable: 

R[e] = md . (2.6) 

One may notice that ^ is a natural variable in the problem since, with A^ = 1, it gives 
the deviation from the critical point: 

With the new variable and 

F{e) = (e^ - 1) T(e^/2, r) , <f{d) = r'^e''^ - e""^^) , (2.8) 
the recursion relation (2.5) reads: 

oo 

F(9) = ip{e) + r-^F{m9) = ^ r''^^ ip{m^ d) (2.9) 

i=o 

where the last expression follows from the iteration process. 
Making use of the expansion 

k=i 

after a straightforward calculation [3] reproduced in appendix 1, one obtains F(9) as 
the sum of a regular and a singular part: 

F(^) = F,,g(^) + F,i,g(^), 



(2.11) 



+ 00 



j=-oo 
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where: 

oo , 2 

^re.n(^) = E^fc^'' ^ ^ " (2.12) 
A;=n 

The singular part, which has the form given in (1.1) with a = and ji = m, 
behaves as: 

^sing(^) = MO) ^^ ^ = • (2-13) 

where, according to equations (2.4-5) and (2.7-8), the exponent x is related to the 
surface magnetization exponent through [16] 

= ^ = - - ^ (2.14) 
2 2 Inm ^ ^ 

when -Fging(^) provides the leading contribution to F{9), i. e. when x < 1 (see below). 

The amplitude A{$), which is log-periodic in 9 with period Inm, may be written 
as the Fourier expansion 



m = J2 M exp (27ri. ^) = 0- Yl r-^'^Prernim^O) . (2.15) 



+ 00 



Inm/ 

j=-oo 



The last relation can be inverted and after some algebra, the Fourier coefficients are 
obtained as [3]: 



+ 00 

27ri3 



M = ^ I du u-'—f^^,^M • (2.16) 







The system displays two different regimes depending on the value of the 
modulation amplitude r. 

When r > r^ = ^Jra^ from (2.13) one obtains x>\ and the regular part in (2.11) 
gives a leading contribution to Fid) which is linear in d. This linear term corresponds 
to a ^-independent leading contribution to T(A, r) according to (2.8) and 



5(A,r) = 2+^ ^ A{Q) Q""-^ ^ 0{Q) . (2.17) 

— m 



It follows that there is surface order at the critical point with m^^ ~ (r — r^^l'^ and 
a first-order surface transition [16]. The deviation from the critical magnetization 
behaves as d^~^ with an oscillating amplitude for ^Jm < r < m and is linear in 6 
above. 

When r takes the form m"^/^ so that x is equal to the integer n, a logarithmic 
behaviour is obtained. This is the case in particular at when the surface transition 
changes from first to second order. Let us write x = n — e and extract from F^^^{9) 
the term k = n which is now singular, so that: 

m = Ke,{0) - ^ + A{e) r (1 - 6ln^) . (2.18) 



Figure 1. Variation of the leading amplitude Aq with the modulation amplitude r 
for the Fredholm sequence with m = 2. It diverges at Vc = when the surface 
transition changes from second to first order (see (2.19)). 



The leading contribution to the amphtude A{9) comes from Aq, the first term in 
V'rem(^) leading to a singular contribution at the lower hmit in (2.16). Introducing a 
cut-ofi^ at ^Q, one obtains: 

4 ^ ^ ['°du u-'+^ = p^^ (2.19) 
In m 7o In m e 

which finally gives: 

m = Ke^io) - r + o(r) . (2.20) 

^ mm 
As a consequence, at the magnetization vanishes as 

- lnm\ 



= (^) ' (-ln9)-V^+0(l). (2.21) 
V m — 1 / 



When 1 < r < r^, which corresponds to < a; < 1, the remainder starts on n = 1 
and since (/Jq in (2-10) vanishes, we have V'rem(^) = fi!^)- Inserting (p{9) given by (2.8) 
in the amplitude in (2.16), after an integration by parts, one obtains 



(e""^" - e"^) u 
A^ = ' 



In + 27ris 



+ 00 _ I 

+ / du ^ ^ 



/" 

^0 



... Inr^ + 27ris 

(2.22) 

1 — /'^^ ^/ Iny^- — 27ris\ 

In + 27ris V In m / 

The oscillating amplitudes involve the Euler gamma function of a complex argument 
T{z). In this regime, as well as for smaller values of r, the behaviour of the surface 
magnetization is governed by F^i^^iO). It gives a contribution 9^~^ to 5" which 
is divergent at the critical point. The surface transition is second order and the 
magnetization exponent is given by (2.14). 



Figure 2. Log-periodic amplitude A(9) as a function of the deviation from the 
critical coupling, = InX^ for the Predholm sequence with m = 2. The oscillation 
amplitude is quickly decreasing when r increases (r = .5, .6, .7, .8). 



Finally when r<l,a;<0, ?i = and once more <^i.em(^) = fi^)- The leading 
power 9^ has an oscillating amplitude and one obtains: 



+00 



-ILL r 

In m /n , ^ 

'^^ (2.23) 

— 1 ^ /lnr~^ — 27ris\ 
Inm \ Inm / 

Explicit expressions for the oscillating amplitudes are given in appendix 2. The 
behaviour of Aq as a function of r for the Fredholm sequence with m = 2 is shown in 
figure 1 whereas the log-periodic oscillations are shown in figure 2 for different values 
of the modulation amplitude. 



3. Critical surface magnetization with a relevant aperiodic sequence 

We now consider a sequence generated through substitution on the digits 1 and with: 

m 

1 ^ s{\) =TT^^~\ ••• 

(3.1) 

^ 5(0) = p • • • • • • OOP 

V 

Starting on 1 with p = 3 and m = 2, the substitution process generates the following 
sequence: 

IIOIIOOOOIIOIIOOOOOO--- (3.2) 

The /;,s satisfy /p^^./ = fk+\ for I = l,m and vanish for I = m + l,p which leads to: 

{mrii^ + l fy.j^^ for I — l,m — 1 
(3-3) 
m tij^j^^ for L = m,p 



Figure 3. Surface magnetization rris as a function of A^^ for tlie relevant aperiodic 
sequence with p = 4, m = 2, and r = 1.4, 1.2, 1, .8, .6, .4 from top to bottom. 
The surface transition is first-order when r > 1 and continuous with an essential 
singularity when r < 1. 



The last relation with I = p can be iterated to give = on a sequence with length 
L = p' . The density of defects n^^/L = {m/py vanishes asymptotically as L^~^ where 

is the wandering exponent of the sequence. It follows that the aperiodic modulation 
does not change the critical point of the Ising quantum chain and its bulk critical 
behaviour. For m > 1, the wandering exponent is positive and, according to the 
Luck's criterion, the sequence gives a relevant perturbation. As a consequence one 
may expect a change in the surface critical behaviour. 

Some general results are known about the critical behaviour of the surface 
magnetization in the case of a relevant aperiodic perturbation [23] . When the couplings 
are weakened near to the surface (here for r < 1), as shown in figure 3 the surface 
magnetization vanishes with an essential singularity: ~ exp[— c where t 

is the deviation from the critical coupling and uj the wandering exponent given in (3.4). 
When they are strengthened (r > 1) there is surface order at the bulk critical point 
and the critical magnetization vanishes at = 1 like 

msc-{r-l)^- (3.5) 



Let us now calculate the oscillating amplitude of the critical surface magnetization^ 
by considering the recursion relation for T(r) which is defined as 



oo 



T{r) = 5(A^, r) - 1 = ^ r'^"^' . (3.6) 

Writing — the sum over j is replaced by a double sum: one over /c = 0, oo and 

the other over I = l,p. Making use of (3.3) one obtains the one-parameter recursion 

nr) = ^^-:^^+pT{n. (3.7) 

t Outside the critical point, a two-parameter recursion relation would be obtained. 



Figure 4. Amplitude of the critical magnetization function of the deviation 

6 = Inr^ from rc = 1 for the relevant apriodic sequence with p = 4 and m = 2. 
Near Vc one obtains log-periodic oscillations around Aq ^^"^ = Vl21n2/7r. The 
influence of the correction terms in (3.12) are visible far from Tc- 



With ^ = Inr^ ~ r — r^, we recover the hnear renormalization (2.6). The recursion 
relation (3.7) may be rewritten as (2.9) with replaced by p and 

^ oo 

m = T(e^/^) , m = ^,- =Y.^uO'. (3.8) 

fe=0 

The remainder V'reml^)' which contains the terms in (p(9) such that pm'^ > 1, is <f{9) 
itself so that we have: 

k=o ^ 

(3.9) 

j=-oo 



The singular part diverges as A{9) 9^ with x — — Inp/lnm and the critical 
magnetization behaves as: 

TTig^ {r — r^^"^^ (3.10) 

in agreement with (3.4-5). 

Finally, with (^rem(^) = fi'^) given by (3.8), the Fourier coefficients in (2.16) are 
now given by: 



1 /■ + °° lnp-27ri3 /I m \ 

im Jo ^e" - 1 e^"^ - 1/ 



'-I-00 

mm 



Inm 

fe=i 



^ / d« ^-1+^^^^^ (e-^^- - me-'^--) (3.11) 







p — m ^ /Inp — 27ris\ ^f^^P ~ 27ris\ 
plnm \ Inm / \ Inm / 
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where C{z) is the Riemann zeta function. The complete expression of the log-periodic 
amplitude is given in appendix 2. 

Collecting these results, one obtains 

r) = A{e) e-^ + ^^^^ + 0{e) . (3.12) 

The amplitude of the surface magnetization showing the influence of the correction 
terms is given in figure 4. 



4. Conclusion 

Some exact results have been obtained for the oscillating amplitude of the surface 
magnetization of two quantum Ising chains with an aperiodic modulation of the 
couplings. The two sequences lead to surface extended perturbations which do not 
change the bulk critical behaviour. 

For the marginal generalized Fredholm sequence, there is surface order at the 
critical point for r > r^ = ^frn. The critical magnetization vanishes as (r — r^)^/^. The 
deviation from the critical magnetization, behaving as t^'^ with = 2(lnr/lnm) — 1, 
has an oscillating amplitude. At r^, the magnetization vanishes as (— Int)"-*^/^. 
When r < r^, the transition is second order with a continuously varying exponent 
(5g = (1/2) — (Inr/lnm) and a log-periodic amplitude involving the gamma function 
of a complex argument. 

With the relevant sequence, the surface is ordered at the critical point for r > 
r^, = 1. The critical magnetization vanishes as (r — r^)^/^^"^) where lo = Inm/lnp 
is the wandering exponent of the aperiodic sequence. The amplitude of the critical 
magnetization is log-periodic. It contains a product of gamma and zeta functions of 
the same complex argument. For r < 1, the surface magnetization vanishes with an 
essential singularity as a function of t. 

These aperiodic quantum chains may be considered as discrete realizations of the 
Hilhorst-van Leeuwen model [24] for which the couplings, A^, = A(l -\- g/K^), decay 
continuously as a power of the distance from the surface. The decay exponent y 
corresponds to 1 — w, where u is the wandering exponent of the sequence. 

For both sequences the leading surface critical behaviour is the same as for the 
Hilhorst-van Leeuwen model [22,24] (with the correspondence g — > Inr/lnm for the 
marginal sequence [16]). The difference lies in the occurence of log-periodic critical 
amplitudes for the aperiodic systems, which is linked to the explicit breaking of the 
continuous dilatation symmetry introduced by the modulations. 

Finally, one may notice that for the Fredholm sequence, the critical magnetization 
does not show any oscillating amplitude. The same is true for the Rudin-Shapiro 
sequence which leads to a relevant perturbation [23]. Although discrete scale 
invariance is necessary for the occurence of log-periodic critical amplitudes, it is not 
sufficient. 
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Appendix 1. Splitting F{9) into regular and singular parts 

Inserting the series expansion (2.10) for (p{9) into (2.9) gives: 

oo oo 

F{e) = J2r-''J2^km^'e'' (Al.l) 

j=0 k=l 

The order of the two sums can be changed when the sum over j, which is a geometric 
series, converges. This occurs for r~'^m^ < 1, i. e. for A; < n defined in (2.12). 
Rewriting (p{d) as: 

n-l 

m-Y.^kO^ + ^.er.{0), (A1.2) 

one may spht F{9) in two parts such that 



fc=i " ^ " j=o 



= E 1 lU^U + E-"'Vren.(m^'^) (Al.3) 



where the remainder, v'rem(^) given in (2.12), contains that part of the expansion 
leading to a divergent geometric series in (Al.l). The last sum can be extended from 
— oo to +00 by adding and substracting 

— 1 +00 oo oo /j/j 

E .-v„„(m^e) = E E = - E 137^? ■ f^^-^) 

j= — oo j—1 k=n k=n 

In this way one obtains: 

oo +00 

= E 1 ^'-2^k + E --'^■^ren.K^) , (Al.5) 
k—1 j= — oo 

leading to (2.11). 



Appendix 2. Explicit results for the oscillating amplitudes 



,10 = 



Let us write the Euler gamma function of a complex argument as 

With the Predholm sequence, when 1 < r < equations (2.15) and (2.22) lead 

to: 



A^^s A v-^ (1 — r ^) / ln6' TTs \ 

A(d) = ^„ + > cos 27TS ha- arctan ; — , 

t^, VlnV + 7r2s2 V Inm ' InrJ 



(1 — r ) Ps _ InO ^ _ _^ TTS 

(A2.1) 

21nr 



mm/ 
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with [25]: 



1 2 cxj 

''.-(-l:^)n 



1 + 



«s = 1 V' 1- 

Inm 



fe=o L 



27r,s 



Inm 



fe=0 



(/e + 1) Inm — Inr^ 
27rs 



-1/2 



(A2.2) 



(/c + 1) Inm— In 



— arctani 



271S 



\ (/c + 1) Inm— Inr^ 



where i/' is the digamma function. 

When r < 1 equations (2.15) and (2.23) give: 



A{9) = Aq + 2 



r-^-1 
Inm 



27TS 



s=l 



In^ \ 

] + «s ) ' 

mm / 



where now: 



r ^ — l_/lnr ^ 
° Inm \ Inm 



" Vli^^J 11 ^ Ulnm + lnr-2j 

27rs 



(A2.3) 



27rs , /Inr 

<^s = ~i 

mm \ mm 



k=0 

2 , °o 



+E 



/c In m + In r ^ 



-1/2 



arctan 



27rs 



(A2.4) 



/c In m + In r ^ 



With the relevant sequence, the Fourier expansion in (2.15) together with (3.11) 
gives: 



oo oo 



plnm \lnm/ \lnm/ 



(A2.5) 



where: 



Xsk = 27rs 



In/c 
Inm 



(A2.6) 



whereas Pg and are now given by (A2.4) with r replaced by p. 
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